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Abstract 

The first observation of about 2000 candidates, with a background contamination below 3%, 
of the rare decay K ± —> 7 r ± 7 r°e + e _ is reported by the NA48/2 experiment. The preliminary 
branching ratio in the full kinematic region is obtained to be: —> TT ± n°e + e~) = (4.06 ± 

0.17) • 10 -6 by analyzing the data collected in 2003. A sample of 4.687 x 10 6 
decay candidates with a negligible background contamination collected in 2003-04 is analyzed 
to search for the dark photon ( A') via the decay chain K ± —» 7 r ± 7 r°, 7 r° —> yA', A' —> e + e~. No 
signal is observed, and preliminary limits in the plane dark photon mixing parameter e 2 versus 
its mass tua 1 are reported. 


1 The NA48/2 experiment 

The NA48/2 experiment at the CERN SPS collected a large sample of charged kaon {K±) decays 
during its 2003-04 data taking period. The NA48/2 beam line has been designed to deliver simul¬ 
taneous narrow momentum band K + and K~ beams originating from the collision of the primary 
400 GeV/c protons extracted from the CERN SPS on a beryllium target. Secondary beams with 
central momenta of (60 ± 3) GeV/c (r.m.s.) following a common beam axis were used. The beam 
kaons decayed in a fiducial decay volume contained in a 114 m long cylindrical vacuum tank. The 
momenta of charged decay products were measured in a magnetic spectrometer, housed in a tank 
filled with helium placed after the decay volume. The spectrometer comprised four drift chambers 
(DCHs) and a dipole magnet. A plastic scintillator hodoscope (CHOD) producing fast trigger 
signals and providing precise time measurements of charged particles was placed after the spec¬ 
trometer. Further downstream was a liquid krypton electromagnetic calorimeter (LKr), an almost 
homogeneous ionization chamber with an active volume of 7 m 3 of liquid krypton, 27A'o deep, seg¬ 
mented transversally into 13248 projective ~2x2 cm 2 cells and with no longitudinal segmentation. 
An iron/scintillator hadronic calorimeter and muon detectors were located further downstream. A 
dedicated two-level trigger was used to collect three track decays with a very high efficiency. A 
detailed description of the detector can be found in [Tj. 

2 First observation of K ± -A 7r ± 7r°e + e _ decay 

The K ± —> 7r ± 7r°e + e decay proceeds through virtual photon exchange which undergoes internal 
conversion into electron-positron pair, i.e. K± —> 7 r =t: vr 0 7 *—> 7 r = 1 : 7 r 0 e + eA The 7 * is produced 
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by two different mechanisms: Inner Bremsstrahlung (IB), where the 7 * is emitted by one of the 
charged mesons in the initial or final state and Direct Emission (DE) when 7 * is radiated off at the 
weak vertex of the intermediate state. As a consequence the differential decay width consists of 
three terms: the dominant long-distance IB contribution (pure electric part E), the DE component 
(electric E and magnetic M parts) and the interference between them[2j. The interference term 
collects the different contributions, IBE, IBM and EM. For this reason the 7 r =t 7 r°e + e~decay offers 
interesting short and long distance parity violating observables. In the K± —> 77 ^ 77 ° y mode the 
interference consists only of the IBE term[3j, because the remaining (EM) interferences are P- 
violating, but cancel out upon angular integration. There are few theoretical publications related 
to the 77 ^- 77 ® e + e-mmm- Recently authors of [2j where able to predict, on the basis of the 

NA48/2 measurement of the magnetic and electric terms in K± —>• 7 r ± 7 r 0 7 [ 6 ], the branching ratio 
of the single components. No experimental observation has so far been reported. 

2.1 Selection and background estimates 

A± —> 7r ± 7r°e + e event candidates are reconstructed from three charged tracks and two photons, 
forming neutral pion, pointing to a common vertex in the fiducial decay volume. Particle identifica¬ 
tion is based on the energy deposition in LKr (E) associated or not to a charged track momentum 
(p) measured in the spectrometer. The charged track is identified as electron/positron if its E/p 
ratio is greater than 0.85, and as a charged pion if the E/p ratio is lower than 0.85. Two isolated en¬ 
ergy clusters without associated track in the LKr are identified as the two candidates photons from 
the 7 T° decay. Their invariant mass is required to be within ±10 MeV/c 2 from the nominal PDGjT] 
7 T° mass. The reconstructed invariant mass of the 7 r ± 7 r°e + e _ system is required to be within ±10 
MeV/c 2 from the nominal PDG[7j K± mass. Two main sources of background are contribution to 
the signal final state: K± 77±77 ®t7® d (A^d) when one of the photon is lost, and —> 77 ± 77d{i) 
(K 2 w d), where 77 Q d denotes the vr° Dalitz decay vr° —> e + e~'y. The suppression of the K^d back¬ 
ground events is obtained by requiring the squared invariant mass of the n + n 0 system to be greater 
than 120 MeV 2 /c , exploiting the presence of three particles with almost the same mass in the final 
state. In order to reject A?rZ) background contamination both the invariant masses M eeil 2 are 
required to be more than 7 MeV/c 2 away from the nominal mass of the neutral pion. Analyzing the 
2003 data, a sample of 1916 signal candidates has been selected with a background contamination 
below 3%. In particular MG simulation predicts a contribution of (26±5.1) candidates form AAd 
and (30±5.5)from AwD events. The normalization mode (A^d) is recorded concurrently with 
the signal mode , using the same trigger logic. A common event reconstruction is considered as 
much as possible aiming to cancel of systematic effects such as particle identification and trigger 
inefficiencies. The selection of the normalization mode AttZ? uses the same set of requirements as 
the signal selection except for the 7 r°-reconstruction and background suppression parts. The neutral 
pion is reconstructed by requiring only one 7 -candidate cluster and computing its invariant mass 
with the electron and positron pair. The only background source for the normalization channel is 
the 3 D mode ( K ± -» ji + v 77 Q d ). In the whole 2003 data sample 6.715 million AAz? candidates 
are selected with a background contamination smaller than 0 . 1 %. 

2.2 Branching ratio measurement 

The total Branching Ratio of K± —> 7 r ± 7 r°e + e _ is obtained using the expression: 

B(K± -> rVe+e-) = ~ Nb Y" B(N) (1) 

Ahv A s es 


2 



where Ns,b,n are the number of signal (1916), background (55.8T7.4) and K^jj events. As,n £S,n 
are the acceptances and trigger efficiencies of the signal and normalization modes. The normal¬ 
ization mode branching ratio B(N) = (2.425 ± 0.076) • 10~ 3 is obtained from the PDGfTj world 
average. The trigger efficiencies (e), very similar (~98%) for signal and normalization mode, are 
measured on data using control samples. The acceptances of the signal, the normalization and the 
background channels are computed using GEANT3-based|8] MC simulations which include the full 
detector and material description, stray magnetic fields, beam line geometry and local detector 
imperfections. The MC simulation for the different K± —> 7r ± 7r°e + e _ contributions IB, DE, and the 
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Figure 1: Reconstructed 7r =1: 7r 0 e + e _ invariant 
mass distributions of the data and simulated 
background samples. 
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Figure 2: Reconstructed e + e _ invariant 
mass distributions of the data and simulated 
background samples. 


electric interference, have been generated separately according to the theoretical description given 
in [2] neglecting the magnetic interference in the present preliminary result. Due to limited statistic 
of the data sample the extraction of the DE and electric interference is not possible in this analysis. 
The signal acceptance has been obtained from a weighted average of the single components accep¬ 
tances, using as weights the relative fractions computed in [2] on the basis of the measurement of 
magnetic and electric terms of —> 7r ± 7r°7 in [6]: 
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To take into account the E,M measurement uncertainties (6], the weights entering the total signal 
acceptance were varied accordingly resulting in a ~ 1% relative change quoted as the systematic 
uncertainty due to the acceptance modeling. As radiative corrections to the K± —» 7r =1: vr 0 e + e~mode 
are not computed in [2], the NA48/2 signal MC simulation included the following effects: the 
classical Coulomb attraction/repulsion between charged particles and the real photon(s) emission 
as implemented in the PHOTOS package The preliminary result for the total branching ratio is 
obtained: 

B(K± -> 7^7r°e + e -) = (4.06 ± 0.10 rtat . ± 0.06 S2/St . ± 0.13 ext .) ■ 10“ 6 (3) 

where systematic errors include uncertainties on acceptance, particle identification, trigger effi¬ 
ciencies and radiative corrections. The external error originating from the normalisation mode 
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Figure 3: The K ± — > 7r ± 7r°e + e preliminary branching ratio is plotted with its experimental error 
(shaded blue band) and its total error (shaded green band). 


branching ratio uncertainty is the dominant error in the present measurement obtained with an 
overall precision of about 3%. The comparison with theoretical expectations is presented in Figure 
[3} The small dashed blue line represents the theoretical prediction with no isospin breaking cor¬ 
rection published in [2]. The big dashed blue line shows the expected isospin breaking corrected 
branching ratio (private communication from the authors of [2j). The experimental value of the 
B(K± — > 7r =l: 7r 0 e + e _ ) is in a very good agreement with the theoretical predictions (within one 
standard deviation). The NA48/2 data sample analyzed has no sensitivity to the DE and INT 
contributions to the M ee spectrum within the current statistics (see Figure [ 2 ]). It will be difficult to 
perform a full Dalitz plot analysis without a proper description of the radiative effects, particularly 
relevant in a final state with two electron/positron. 

3 Search for the dark photon in 7T° decays 

The large sample of 7 r° mesons produced and decaying in vacuum collected by NA48/2 allows for 
a high sensitivity search for the dark photon (A 7 ), a new gauge boson introduced in hidden sector 
new physics models with an extra 17(1) gauge symmetry. In a rather general set of models, the 
interaction of the dark photon (DP) with the ordinary matter is through kinetic mixing with the 
Standard Model hypercharge 17(1) [9]. In these models, the new coupling constant e is proportional 
to the electric charge and the dark photon couples in exactly the same way to quarks and leptons. 
These scenarios could provide an explanation to the observed rise in the cosmic-ray positron fraction 
with energy, and could offer an explanation to the muon gyromagnetic ratio (g — 2 ) anomaly [TOJ. 
The simplest DP model is characterised by two free parameters, the DP mass uia 1 and the mixing 
parameter with the standard model e. Its possible production in the 7 r° decay and subsequent 
decay proceed via the following chain: , 7 r° —> 7 A 7 , A 7 — > e + e _ , producing a final 

state with three charged particles and a photon. The expected branching fraction of the 7 r° decay 
is pi]: g 

I3(7r 0 —> 77 ), (4) 

m W 


B(tt° -> 7 A 7 ) = 2 e 2 1 - 
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with a strong kinematic suppression of the decay rate for DP masses approaching m n o . In the mass 
range 2 m e <C 7774 / < ° accessible in this analysis, assuming that the DP can only decay into 

SM fermions, B(A' — > e + e~) ~ 1 while the allowed loop-induced decays ( A' —> 37 , A' —>• vu) are 
highly suppressed. The maximum DP mean path m in the NA48/2 experiment corresponds to 
an energy of approximately E max = 50 GeV: 

. „ /l0- 6 \ /100 MeV\ 2 

Tmax ~ ( E ma , x /m A ')cT ~ 0.4 mm x I - I x ^— ——— j , 

In the accessible parameter range (m^ > 10 MeV/c 2 and e 2 > 5 x 10~ 7 ) L max does not exceed 
10 cm and the DP can be assumed to decay at the production point. In this prompt decay scenario 
the NA48/2 3-track vertex reconstruction does not introduce significant acceptance losses as the 
typical resolution on the vertex longitudinal coordinate is ~ 1 m. The DP signature is identical to 
that of the Dalitz decay —> e + e~ r y , which therefore represents an irreducible background and 
limits the sensitivity. The largest ir^ sample, and therefore the largest sensitivity, is obtained form 
the study of the decays. 

3.1 Event selection and background simulation 

The full NA48/2 data sample is used for the analysis. The K^d event selection requires a three- 
track vertex reconstructed in the fiducial decay region formed of a pion (vr^) candidate track and 
two opposite-sign electron (e^) candidate tracks. Charged particle identification is based on the 
ratio of energy deposition in the LKr calorimeter to the momentum measured by the spectrometer, 
which should be smaller (greater) than 0.85 for pion (electron) candidates. Furthermore, a single 
isolated LKr energy deposition cluster is required as the photon candidate. The reconstructed 
invariant mass of the system (Fig. [4| is required to the consistent with the K± mass. A 

sample of 4.687 x 10 6 fully reconstructed tt g d decay candidates in the e + e~ invariant mass range 
m ee > 10 MeV/c 2 with a negligible background is selected. The candidates mainly originate from 
K-i-kD decays, with 0.15% coming from the semileptonic decays (denoted K^d 

below). Correcting the observed number of candidates for acceptance and trigger efficiency, the 
total number of decays in the 98 m long fiducial decay region for the analyzed data sample 
is found to be Nk = (1.55 ± 0.05) x 10 11 , where the quoted error is dominated by the external 
uncertainty on the tt® d decay branching fraction The reconstructed e + e~ invariant mass 

(m ee ) spectrum of the K' 2 -kD candidates is displayed in Fig. [5] A dark photon produced in the 
7 x° D decay and decaying promptly to e + e _ would appear as narrow peak in the spectrum. Monte 
Carlo (MC) simulations of the and K^r> processes are performed to subtract the irreducible 

7 t° d background. The 7 x° D decay is simulated using the lowest-order differential decay rate in 12j. 
Radiative corrections to the differential rate are implemented following the approach of Mikaelian 
and Smith |T2j recently revised to improve the numerical precision H3j. The method introduces 
only weights 5(x,y) and does not account for the emission of inner bremsstrahlung photons. 

3.2 Dark photon search technique 

A search for the DP is performed assuming different mass hypotheses with a variable mass step. 
The mass step of the scan and the width of the signal mass window around the assumed DP mass 
are determined by the resolution on the e + e _ invariant mass. The mass step of the DP scan is 
set to be <r m /2, while the signal region mass window for each DP mass hypothesis is defined as 
±1.5cr m around the assumed mass (both the scan step and the mass window half-width are rounded 
to the nearest multiple of 0.02 MeV/c 2 ). The mass window width has been optimised with MC 
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Figure 4: Reconstructed vr^vr^ invariant 
mass ( 771271 -) distributions of the data and 
simulated background samples. The selec¬ 
tion condition is illustrated with arrows. 



Figure 5: Reconstructed e + e _ invariant 
mass distributions of the data and simulated 
K 27T d and K^ D samples. 


simulations to obtain the highest sensitivity to the DP signal, determined by a trade-off between 
7 t° d background fluctuation and signal acceptance. 

In total, 398 DP mass hypotheses are tested in the range 10 MeV/c 2 < m ee < 125 MeV/c 2 . The 
lower limit of the considered mass range is determined by the limited precision of MC simulation 
of background at low mass, while at the upper limit of the mass range the signal acceptance drops 
to zero. The numbers of observed data events in the signal region (IV 0 bs) and the numbers of 
iTp background events expected from MC simulation corrected by the measured trigger efficiencies 
(IVexp) in the DP signal window for each considered mass hypothesis are presented in Fig. [6j They 
decrease with the DP mass due to the steeply falling -k^ differential decay rate and decreasing 
acceptance, even though the mass window width increases, being approximately proportional to 
the mass. 

Confidence intervals at 90% CL for the number of A' —> e + e~ decay candidates (N^p) in each 
mass hypothesis (IVdp) are available from !V 0 b s , N exp and 5N exp using the Rolke-Lopez method |14j 
assuming Poissonian (Gaussian) errors on the numbers of observed (expected) events. For the 
preliminary results, it is assumed conservatively that !V 0 bs = N exp in cases when -/V 0 b s < N exp , 
as the employed implementation of the method (from the ROOT package) has been found to 
underestimate the upper limits in that case. Upper limits at 90% CL on B(n° —> 7 A') in each DP 
mass hypothesis in the assumption B{A' —> e + e - ) = 1 are computed using the relation 


R(vr° -7 jA') = \B(K 2n )A(K 2n ) + B(K^)A{K^ 


N k L 


The acceptances A(K 2n ) and A{K ^ 3 ) of the employed K 27T d event selection for the I\ 2n and 
decays, respectively, followed by the prompt 7r° — > 7 A', A' —> e + e _ decay chain, are evaluated 
for each considered DP mass with MC simulation. Event distributions in the angle between e + 
momentum in the e + e~ rest frame and the e + e _ momentum in the 7T° rest frame are identical for 
the decay chain involving the DP (7r° —> ^A', A' — > e + e ~) and the decay, up to a negligible 
effect of the radiative corrections that should not be applied in the former case. Therefore DP 
acceptances are evaluated using the MC samples produced for background description. 

The largest uncertainty on the computed B(ir° — > 7 A!) is the external one due to B(ir°) entering 
via Nk- It amounts to 3% in relative terms and is neglected. The obtained upper limits on 
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Figure 6: Numbers of observed (red) and ex¬ 
pected (blue) DP candidates in the DP signal 
mass window as for each assumed DP mass 
hypothesis. 


Figure 7: The NA48/2 preliminary upper 
limits at 90% CL on the mixing parameter 
e 2 versus the DP mass compared to the 
other published exclusion limits. 


B(tt° —> 'yA') are 0 ( 10~ 6 ) and do not exhibit a strong dependence on the assumed DP mass, as the 
negative trends in background fluctuation (Fig. [6]) and acceptance largely cancel out. Upper limits 
at 90% CL on the mixing parameter e 2 in each considered DP mass hypothesis are calculated from 
those on £>( 7 r° —> 'yA') using Eq. (jd]). The resulting preliminary DP exclusion limits, along with 
constraints from other experiments m, the band of phase space where the discrepancy between 
the measured and calculated muon g — 2 values falls into the ±2 a range [TUI Hi] due to the DP 
contribution, and the region excluded by the electron g — 2 measurement, are presented in Fig. [7j 
The obtained upper limits on e 2 represent an improvement over the existing data in the DP mass 
range 10-60 MeV/c 2 . Under the assumption that the DP couples to SM through kinetic mixing 
and decays predominantly to SM particles, the NA48/2 preliminary result excludes the DP as an 
explanation for the muon (g-2) anomaly in the range 10-100 MeV/c 2 . 
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